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It has been pointed out that it is impossible to obtain a unitary implementation of the dynamics for the 
polarized Gowdy T 3 cosmologies in an otherwise satisfactory, nonperturbative canonical quantization proposed 
for these spacetimes. By introducing suitable techniques to deal with deparametrized models in cosmology 
that possess an explicit time dependence (as it is the case for the toroidal Gowdy model), we present in this 
paper a detailed analysis about the roots of this failure of unitarity. We investigate the impediments to a unitary 
implementation of the evolution by considering modifications to the dynamics. These modifications may be 
regarded as perturbations. We show in a precise manner why and where unitary implementability fails in our 
system, and prove that the obstructions are extremely sensitive to modifications in the Hamiltonian that dictates 
the time evolution of the symmetry-reduced model. We are able to characterize to a certain extent how far the 
model is from unitarity. Moreover, we demonstrate that the dynamics can actually be approximated as much as 
one wants by means of unitary transformations. 
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I. INTRODUCTION 

Symmetry-reduced models have been used over the past 30 
years as an appropriate arena to test strategies aimed for the 
quantization of the full theory of gravity within the canon- 
ical approach, as well as toy models which can provide us 
with insights about the kind of phenomena that should be ex- 
pected in quantum general relativity. Most of the examples 
of symmetry-reduced models studied so far are minisuper- 
space models 1 1 ], namely, simple systems where the reduction 
leaves only a finite number of physical degrees of freedom. 
There is another class of models which has more interest inas- 
much as they retain the field complexity of general relativity. 
These are the so-called midisuperspace models (for a recent 
review see Ref. |2]), which after reduction possess an infinite 
number of degrees of freedom. Thus, their quantization would 
lead to a true quantum field theory. 

Within this class, and together with the Einstein-Rosen 
waves |3|, the model that has deserved more attention lately 
is the Gowdy r 3 cosmological model | % |2, |y, la, la] . This 
model was introduced by Gowdy during the seventies in a sys- 
tematic search for all spacetimes with two commuting space- 
like Killing vector fields and compact spatial hypersurfaces 
11(11 . Apart from the T } topology of a three-torus, the other 
two possible spatial topologies for the Gowdy spacetimes are 
the three-handle S 1 x S 2 and the three-sphere S 3 (or the topol- 
ogy of a manifold covered by one of the above). The interest 
in the Gowdy T 3 model can thus be easily understood, since 
it provides the simplest of all the inhomogeneous, empty, 
spatially closed cosmological systems. The genuine field- 
theory character of this model and its possible applications in 
cosmology make it a natural candidate to study fundamental 
questions about canonical quantum gravity and quantum field 
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theory in curved spacetimes. 

Its quantization was already considered in the seventies 
fllll . and revisited both in the eighties [H, Q and 
more recently 0, 0, • The first preliminary attempts 
to define a quantum theory and extract physics from the 
Gowdy r 3 model fill O were followed by more detailed 
analysis that discussed the nonperturbative quantiza- 

tion of the system employing the Ashtekar formulation of 
Lorentzian general relativity Hal . Considerable progress has 
been achieved lately in defining a complete quantization of 
the (sub-)model with linear polarization, in which both Killing 
vectors are hypersurface orthogonal |5]. The proposed quan- 
tization is based on the fact that the polarized model can be 
equivalently treated as 2 + 1 gravity coupled to a massless 
scalar field, defined on a manifold whose topology is T 2 x EL 

One important aspect in the study of quantum cosmolog- 
ical models is their dynamical evolution. For the polarized 
Gowdy T 3 model with the particular quantization performed 
in Ref. |5|], it has been recently shown that the dynamics is 
not implementable at the quantum level as a unitary transfor- 
mation 

El [J. From the point of view of canonical quantum 
gravity, this result does not represent a serious drawback for 
the simple reason that (owing to the compact nature of the 
spatial slices) time evolution is pure gauge in the Hamiltonian 
description. Hence, there is no time evolution and no dynam- 
ics. The system is endowed with a fictitious dynamics via a 
"deparametrization" procedure, and there is no apparent rea- 
son to select a preferred deparametrization. 

Nevertheless, if one accepts that unitary evolution is a key 
ingredient in conventional (field) quantum theory, necessary 
in order to pose physically meaningful questions for issues 
like those concerning the initial singularity in cosmology, the 
lack of a unitary time evolution is a drawback for the kind 
of quantization put forward in Ref. |5]. From that quantum 
theory, one would not be able to extract predictions for differ- 
ent instants of time, because probability is not conserved. In 
this sense, the quantization is not fully consistent |6|. Thus, 
restoration of unitarity in the evolution seems a fundamental 
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issue in order to achieve a satisfactory quantization of the po- 
larized Gowdy r 3 model. A rigorous quantization along these 
lines will provide us with a specific example of midisuper- 
space that can be very helpful, as a point of reference for com- 
parisons, for a future quantization starting from loop quantum 
gravity (where impressive progress has been made, but exclu- 
sively for minisuperspaces 1 16]) or for implementations of the 
"consistent discretization" approach] 1711 (which has recently 
been applied to the Gowdy model 1 lSTI). 

In this work, we will explore the reasons behind the failure 
in the unitary implementability of the dynamics as a first step 
in proposing solutions to it or introducing alternative descrip- 
tions for the quantum evolution. Although the lack of unitar- 
ity certainly follows from the absence of square summability 
for the antilinear part of the Bogoliubov transformation which 
relates the annihilation (creation) operator at, say, the "final" 
time tf with the annihilation and creation operators at an "ini- 
tial" time ?, |6j|7|], we want to analyze in detail the roots of this 
failure. We will show in a precise manner why unitary imple- 
mentability fails in our system and, in a certain sense, we will 
be able to characterize how far the model is from unitarity. In 
doing so, we will introduce suitable techniques to deal with 
deparametrized models in cosmology that possess an explicit 
time dependence, like the Gowdy one. 

The structure of the paper is as follows. In Sec. |n]we ob- 
tain the reduced phase space of the system by performing the 
symmetry reduction and introducing a deparametrization pro- 
cedure in order to (partially) fix the gauge freedom 1 and select 
a Hamiltonian vector field to represent the dynamics. This re- 
duction and gauge-fixing process had not been presented be- 
fore starting with the spacetime metric of the model in general 
relativity, although the resulting description of the vacuum 
Gowdy spacetimes is essentially the same that was discussed 
in Refs. o|6J (except for the remark on footnote 1 ) and Ref. 
|0]. Therefore, the quantum theory on which we will base our 
analysis is that constructed by Pierri |5] (or, equivalently, that 
analyzed in Ref. [7]). 

Sec. |lll]is divided in three parts. In the first one, we review 
the dynamics of the scalar field that represents (most of) the 
true degrees of freedom of the theory. In this part, a crucial 
remark is that the coordinates of the covariant phase space 
(namely the coefficients that determine the field in terms of 
an orthonormal basis of solutions -defined essentially by the 
negative of the complex structure given in Ref. |3]-) do not 
evolve in time. In these coordinates, the generator of the evo- 
lution is obviously the zero Hamiltonian. The dynamics will 
be introduced through a time-dependent map from the covari- 
ant to the canonical phase space. In the remaining parts of the 



The "q-number" nature of the time variable which occurs in Refs. 01 
will be avoided here by introducing a gauge condition that is slightly differ- 
ent to the one imposed in those references; in this way the physical degrees 
of freedom will be neatly disentangled from the time variable. This type of 
mixing is also absent in Ref. 7 ] ; however, the disentanglement is attained 
there thanks to the introduction of an appropriate (partially) reduced line 
element since the very beginning, rather than to its construction by gauge 
fixing. 



section, we will investigate the impediments for a unitary im- 
plementation of the evolution by considering modifications to 
the dynamics that may be regarded as small perturbations. We 
will be able to identify where the failure of unitarity comes 
from and, in the last subsection, prove that the severity of the 
problem is greatly ameliorated by the fact that small correc- 
tions to the dynamics can be implemented in a unitary way. 
In addition, our analysis makes clear that the diagonalization 
of the Hamiltonian performed in Ref. jfj, 0] is just an instan- 
taneous diagonalization which ignores the change in time of 
the Bogoliubov coefficients. Here, these time variations are 
explicit and rigorously taken into account. Finally, the con- 
clusions and some further comments are presented in Sec. IIVI 
One appendix is added which contains a proof about the be- 
havior of the coefficients employed in the main discussion. In 
the following, lower case Latin indices on a tensor will denote 
its purely spatial components, whereas capital case Latin in- 
dices will be used to denote the tensor itself (abstract index 
notation). 

II. THE POLARIZED GOWDY MODEL 

The polarized Gowdy T 3 model describes globally hyper- 
bolic four-dimensional vacuum spacetimes, (AL gAs), with 
two commuting hypersurface orthogonal spacelike Killing 
fields and compact spacelike hypersurfaces homeomorphic to 
a three-torus. Since global hyperbolicity implies that we can 
foliate {M,gAB) by Cauchy surfaces, parametrized by a 
global time function f, then a 3 + 1 decomposition is available 
and the line element can be written 

ds 2 = -N 2 dt 2 + hijidx' + N'dtXdxJ + N j dt) , (1) 

where we choose t e R + , the coordinates in the sections of 
constant time are \x'} := \x l = 8, x 2 = cr, x 3 = 6} with x' e S , 
N and \N'} are, respectively, the lapse function and the com- 
ponents of the shift vector N A , and \hij] are the components 
of the induced spatial metric Iiab- 

In addition, we will impose that (d/dx"^ (a = 2,3) are 
the two spacelike Killing vector fields. Thus, the metric must 
be independent of the coordinates x" . Moreover, perform- 
ing a (partial) gauge fixing along the lines explained in Ref. 
IU9I1 (for pure gravitational plane waves) and Ref. I20ll (for 
cylindrical spacetimes), remembering that the metric func- 
tions must be periodic in 6, and using that {d/dx a ) A are com- 
muting hypersurface orthogonal vector fields, one gets a line 
element for the reduced model of the form 

3 

ds 2 = -N 2 dt 2 + h ge [de + N e dt] 2 + ^ h aa (dx a ) 2 . (2) 

Let us consider now the following change of metric vari- 
ables {hij} h-> {Q a } := {(//, y, r}, defined by 

hee = e 7 ^ , /w = e^r 2 , h ss = . (3) 

Since this change is just a point transformation, the momenta 
P a canonically conjugate to Q" are 
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where p'j - yfh(K'j - Kh'i) are the momenta canonically con- 
jugate to hip we have set the Newton constant G equal to n/4, 
h is the determinant of the induced metric, and {K'i} are the 
components of its extrinsic curvature, with trace equal to K. 
Substituting in Eq. (0 our new set of variables and introduc- 
ing the densitized lapse factor N := N/ *\fh, we arrive at the 
line element 



ds 1 = ^^{-T^dt 1 + [de + N e dt\ 2 ^ 
+ e-f^do 2 +e 2 "'d6 2 y 
whose Einstein-Hilbert action is then given by 



(5) 



= £ dt § 



P T f + P y y + P^ - (NC + N°C e ) 



(6) 



The presence of the remaining first class constraints C and 
Cg reflects the fact that the gauge has been only partially 
fixed. These (densitized) Hamiltonian constraint and momen- 
tum constraint are, respectively, 



C := ~(Pj -2tP t P 7 )+ T -{4t" -2y'r' +Tip' 2 ), 
Cg := P t t' + P y7 ' + P^' - 2P' y . 



(7) 



The prime denotes the derivative with respect to 6. 

Note that, since the spatial slices are compact, there exist 
no boundary contributions to the Hamiltonian. Therefore, the 
total Hamiltonian vanishes on the constraint surface and there 
is no distinction between gauge and dynamics. It is then nec- 
essary to carry out a deparametrization in order to introduce 
dynamics. This deparametrization is accomplished as part of 
a gauge fixing of the model: one imposes suitable conditions 
which, together with the constraints |7}, form a set of sec- 
ond class constraints, allowing the reduction of the system. 
Explicitly, we demand the following gauge-fixing conditions 
(which are a slight modification of those introduced in Refs. 
ii) 



Py + P = , 



c?2 



tp = . 



(8) 



The first of these conditions requires the momentum canon- 
ically conjugate to y to be homogeneous (independent of 0). 
Furthermore, this homogenous part is a constant (of motion) 
p. In this sense, it is worth pointing out that the Poisson brack- 
ets of j>P y l(2n) (i.e. -p) with all the first class constraints Q 
vanish weakly, so that it is indeed a Dirac observable. On 
the other hand, the second of our conditions fixes the met- 
ric function t equal to the global time function t except for a 
rescaling that is constant on shell, though can vary on different 
solutions. Modulo constraints and gauge-fixing conditions, a 
straightforward calculation shows then that 



dOGCt 



- P G' . 



g2 



dOFC^ 



tp 2 F 



(9) 



where the smearing functions F and G on S 1 are, respectively, 
a density of weight -1 and a scalar. Therefore, if F and G' 
are different from zero, these Poisson brackets do not vanish 
provided that p + 0. Thus, we have to restrict all considera- 
tions to the sector of solutions with nonzero p in order to get 
a well-posed fixation. 

The next step in this procedure consists in demanding the 
compatibility of the gauge-fixing conditions with dynamics: 
the total time derivative of g\ and g2 must vanish for some 
choice of N and N e . This derivative is the sum of the Pois- 
son bracket with the total Hamiltonian j> Jif and the partial 
derivative with respect to the explicit f-dependence. Modulo 
constraints and gauge-fixing conditions, we have 



gi = -p (n 6 )' , g 2 = -p + tp 2 N 



(10) 



The requirements gi = and g2 — are then satisfied if N e 
is any function of t and = {tpY 1 . It is worth noticing that, 
while the densitized lapse function is completely determined 
in this process, the shift function is not fully fixed. There 
remains some diffeomorphism gauge freedom, generated by 
the homogenous part of the constraint Cg (after reduction). 
Besides, note that we have to further restrict p to be positive 
in order to ensure the positivity of the lapse function 2 . 

In order to extract the true degrees of freedom, one solves 
the set of second class constraints \C,Ce,g\,gi\, obtaining 



PPr 



up? 



^\- L +tp^' 



2\tp 
py' = := A. 



(11) 
(12) 



By performing a Fourier expansion in 6 of the functions y and 
A (which is possible given the smoothness of the fields on S 1 ), 
it is not difficult to see that identity \\2l allows us to solve for 
all modes of y but the zero mode. More precisely, the Fourier 
coefficients y n are determined in terms of A„ by inpy n = A„. 
Thus, there is still an undetermined coefficient, namely yo, and 
consequently we are left with a global degree of freedom. 

Furthermore, note that integration over S 1 of Eq. H2\ leads 
to the global constraint 

K = -^j>dOP f ij/ =0, (13) 

which is essentially the homogenous part of the constraint Cg. 
Therefore, the diffeomorphism gauge freedom has not been 
entirely removed and the ^-component of the shift vector can- 
not be completely fixed. However, as we have already seen, 
the only allowed dependence of N e is that on t. This type of 
shift can always be absorbed by redefining our angular coor- 
dinate 6 |4]. After our gauge fixing and the absorption of the 



Otherwise, one should consider t 6 R instead of t s R + for the time flow 
vector field to be future directed. 
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shift, the metric becomes 

ds 1 = e^^dt 2 +d6 2 ) + e-' f 't z p 2 dij 2 + e' f 'd6 2 , (14) 

q . V A « eM 



7 = 



where q : = V27ryo is the coordinate canonically conjugate to 
—p (the zero mode of P y / V27r). 

The reduced action for the system (modulo a spurious 
boundary term -pq \' t r ) is 

S r = j"J dt^p[q+ <j)d8tpjj+ <j)dO[Prf + pP T ] 



r 



dtipi] + (hde 



dO — 

2p 



f p2 

^ + tp*' 2 
tp 



I pi 
tp 



(16) 



In the first equality, P T denotes the solution given in Eq. 
(lilt . Thus, S r is a functional on the reduced phase space 
F,., which is coordinatized by (77, p, ip, P^), and where the 
(only nonvanishing) basic Poisson brackets are \p, 77} = 1 and 
{\]/{t, 8), P^(t, 8)) = 6(8 - 8). Note that, owing to the presence 
of the global constraint Jl 3i . the space of physical states does 
not correspond to V r but to a submanifold of it. However, 
since this submanifold is nonlinear, the reduction by the con- 
straint is usually postponed to the quantum theory, where it is 
imposed as an operator condition on quantum states. 
Let us now perform the canonical transformation 



tp = yfpip , P = 
1 



Q 



-rj + 



2p 



VP' 
dOP^iff, 



P- 



(17) 



In terms of this new set of phase space variables the reduced 
action reads 

S r = f ' dt [pQ + j)d9 [p$ - J^]J - PQ\* , (18) 



where the (reduced) Hamiltonian density is 



( P 2 \ 

— + td> L 
t 



(19) 



Thus, our midisuperspace model consists of a phase space f r 
coordinatized by the canonical pairs (Q, P) and {tp, P$), which 
we will call the global and local degrees of freedom, respec- 
tively. Remember that P is strictly positive. To arrive at a 
true canonical pair of real variables, we could always replace 
(Q, P) with (QP, InP). There also remains a global constraint 
on the system (Ao = 0) which restricts the physical states 
to lie in a submanifold of f r . Note that, given the (Q, P)- 
independence of the Hamiltonian density, these "point parti- 
cle" degrees of freedom are constants of motion. Hence a 



nontrivial evolution may only take place in the field sector 
T = {((p, Pip)}. Since the time evolution affects only the local 
degrees of freedom, we will focus on them in our analysis. 

Varying action J 1 81 with respect to and P^ one gets the 
field equations 



P = t(f>, P^ = t<p" 



(20) 



Hence, we only have to consider all smooth solutions to the 
second-order differential equation 



1 . 







(21) 



in order to specify the classical spacetime metric. Using the 
method of separation of variables, it is not difficult to see that 
these solutions, that we will genetically denote by ip, adopt the 
formJl 

<p(f,S) = ^ J] [A n H Q {\n\t)e M + A* n Hl{\n\t)e- M \ 



neZ, ni=0 



(q + po In f) , 



(22) 



where the symbol * denotes complex conjugation, qo and po 
are constants . Ho is the zeroth-order Hankel function of the 
second kind 12111 and, in order to guarantee pointwise con- 
vergence, the sequence of constant coefficients \A n } has to 
decrease faster than the inverse of any polynomial in n as 
n — > ±00. Expression ( 1221 determines the metric dl4l > (with 
iff = <pl yfp) except for the values of (q, p). One can show that 
its Kretchmann scalar blows up at t — 0, so that there is an ini- 
tial singularity and the global time function t must be strictly 
positive. 

In the field sector, the physical phase space can be alterna- 
tively described by the submanifold obtained by imposing the 
constraint Ao in F, or by that submanifold of the space V of 
the smooth solutions (I22> defined by the constraint 



An 



:= Yj n K A n = 0. 



(23) 



neZ, «#0 



In addition notice that, for the field sector, the reduced action 
(II 81 can be viewed as that corresponding to an axi-symmetric, 
massless, free scalar field propagating in the fictitious flat 
background in three dimensions: 

(f) gAB = -(dt) A (dt) B + (d8) A (d8) B + t 2 (dcr) A (dcr) B . (24) 

Thus, we can identify F with the canonical phase space of 
the scalar field in this background, (At - T 2 x R + ,^gAB), 
whereas the space V of smooth solutions can be considered 
as the covariant phase space of such a Klein-Gordon field. 
Namely, <p and P^, are the configuration and momenta on the 
constant-time section X, of the scalar field tp propagating in 
(M,^ gAB)- Besides, since the fictitious background is glob- 
ally hyperbolic, given a smooth Cauchy surface £, there will 
be a natural isomorphism I[L to ] between the linear spaces F 
and V. In this framework, the analysis of the dynamics of the 
polarized Gowdy T 3 model becomes equivalent to the study 
of the time evolution of the free scalar field. In the next sec- 
tion we review this dynamics and discuss the obstructions to 
its unitary quantum implementation. 
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III. DYNAMICS 

It has recently been shown that the dynamical evolution 
generated by the reduced Hamiltonian H r - j> M'r [see Eq. 
(I19H cannot be implemented as a unitary transformation, nei- 
ther on the kinematical Fock space |6] constructed from V 
with the complex structure associated with the field decom- 
position 11221 . nor in the physical Hilbert space of states [7] 
determined by the kernel of the operator version of the con- 
straint d23i . We want to analyze in detail the reasons behind 
this lack of a unitary implementation and discuss how severe 
the problem is, studying whether small corrections (coming 
e.g. from quantum or perturbative modifications) to the dy- 
namics may suffice to restore the unitarity. 

For the space V, we will employ as coordinates the con- 
stants coefficients of the field decomposition J22i . whereas for 
r we will use a different set that absorbs in its (implicit) time 
dependence all the evolution of the field. We will see that the 
dynamics in F is dictated by H r , whereas that in V is frozen, 
because the considered coefficients are constants of motion 3 . 

For the sake of completeness and clarity, let us remember 
some definitions and make a few remarks that will be useful 
in our analysis. Firstly, we recall that given two field decom- 
positions in different orthonormal bases of solutions, namely 

<p = ^ n A n fn(t,e)+A*j;(t,e)wd<p = Y Jn A„ gn (t,e)+A; g *„(t,6), 

their coefficients are related by a Bogoliubov transformation. 
That is, A„ = £,„ a mn A m +p* mn A* m with Yik( a iku) k -PikP) k ) = Sy 
and Xki^ikPjk -PikUjk) = jH. 

Secondly, as in the classical case, the quantum theory for 
the scalar field can be formulated either in a covariant or in a 
canonical approach. In fact, it is generally known that, by en- 
dowing the space of (smooth) solutions V = {(f) to the Klein- 
Gordon equation with a complex structure / compatible with 
the symplectic structure, one can construct in a canonical way 
the Hilbert space of the quantum theory as a symmetric Fock 
space 3? on which the basic observables of the theory are rep- 
resented as annihilation and creation operators ll23ll . On the 
other hand, if the scalar field propagates in a globally hyper- 
bolic spacetime, e.g. M ~ X x R + , given an embedding Z /o of 
2 as a Cauchy surface in M one gets, from the covariant com- 
plex structure J on V, the induced complex structure Jo on the 
canonical phase space T 12411. Once we have /o> we know how 
to construct the Schrodinger representation which is unitarily 
equivalent to the Fock one |24] and how to pass from one rep- 
resentation to the other B25I1 . In particular, we have the analog 
of the one-particle Hilbert space and, by applying the creation 
operator, one can construct the n (functional) particle states. 

In addition, as a consequence of this unitary relation be- 
tween the covariant and canonical approaches, if a symplec- 
tic transformation is unitarily implementable with respect to 



3 Let us emphasize that the coefficients in Eq. 1221 do not display any time 
dependence, not only explicitly, but also implicitly. Accordingly, the total 
Hamiltonian in V indeed vanishes. Properly speaking, the nonunitarity 
proved in Ref. [ 6] is that of the transformation generated by H, on V, which 
in turns can be seen to imply the nonunitary character of the dynamics in 
T, rather than in V. 



the Fock representation, so is it with respect to the (unitarily 
equivalent) Schrodinger representation (and vice versa). Re- 
call that a symplectic transformation on V (F) will be unitarily 
implementable with respect to the Fock (Schrodinger) repre- 
sentation if the antilinear part of its quantum counterpart is 
Hilbert-Schmidt on the one-particle Hilbert space 1 26] . For a 
quantum transformation of the Bogoliubov type, the Hilbert- 
Schmidt condition reduces to £„ m ^6« m | 2 < 00 ■ 

Consider now the Schrodinger representation constructed 
from /o (the complex structure induced by the covariant one, 
J). Given a symplectic transformation T on T, we obtain 
an induced complex structure J' Q = TJqT 1 and, associated 
with it, a new Schrodinger representation. The annihilation 
and creation operators in this new representation are related 
with the annihilation and creation operators of the former one 
through a Bogoliubov transformation. The representations 
corresponding to J' and Jo are then unitarily equivalent if 
the antilinear part of this Bogoliubov transformation is square 
summable. 

In particular, if T = T( tuh) represents the time evolution 
from 1 1 to ?2 on F, then the Bogoliubov transformation relates 
the annihilation and creation operators at the instants ?i and 
?2. With an appropriate choice of coordinates in F, the sym- 
plectic transformation T( tu t 2 ) acts on the elements of this space 
exactly as the Bogoliubov transformation on their quantum 
counterparts. Thus, in order to elucidate whether the sym- 
plectic transformation is unitarily implementable it suffices to 
analyze the square summability of the antilinear part of T( t] ,,). 
Let us emphasize that this procedure is equivalent to determin- 
ing whether the Schrodinger representations constructed from 
the same F (associated with the embedding £,,) but with the 
distinct complex structures J\ and Ji, = T^ uh )JiT^ h) are uni- 
tarily related. 

In the following, we will focus our discussion on analyz- 
ing the failure of unitary implementability of the symplectic 
transformation that determines the dynamics in F, rather than 
examining the complex structures induced by it, since both 
procedures are equivalent, as we have just seen. 



A. Evolution in the canonical and covariant approaches 

We start by expanding in Fourier series our canonical vari- 
ables d> and P*: 



«gz V2tt „ eZ V2tt 



(25) 



where the coefficients <j>„ and P"^ are (implicit) functions of the 
global time coordinate. From the basic Poisson bracket be- 
tween <p and P$, it is easily shown that {<f>„, P"!} = 6'". There- 
fore, we can equivalently consider as our canonical phase 
space that whose coordinates are the set of (complex) canon- 
ical pairs {(<p n , PT)}n€Z- We will call it T^p^. Note that this 
space can be decomposed as the direct sum of P and r () , 
where P is the subspace of vectors with 0o = — and 
Fo is the span of those vectors whose only nonvanishing com- 
ponents are precisely those corresponding to <po and Pz. For 
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all positive integers m e N- {0}, let us now consider the trans- 
formations 



where 



(*m — 



. P$, <P-m, P/') ' * ( a m, a* m , a_ m , a*_ m ) , 
\m\<p m + iPT \m\<p- m - iP'£ 



(26) 



(27) 



One can check that these transformations are canonical, 
so that {a n ,ia*J = 5 nm . Hence, the canonical phase space 
can be alternatively described by the symplectic vector space 
T fl = To © f , where the coordinates for f are the (complex 
conjugate pairs of) annihilation and creation-like variables 
{(a m ,a^,a_ m ,fl* m )} meN _( ). The dynamics on Y a [as well as 
in F and T^p^J is dictated by the reduced Hamiltonian for 
the polarized Gowdy T 3 model, H r = £ J4? r , where Jf° r is the 
reduced Hamiltonian density dl9> . Using expressions ( I25i and 
(I27i . and defining 





._ K) 2 






n 




t 2 + 1 








2t 




f 2 -l 








+ It 



m (a m a* + a_ m a_ m ) 



m (a* a* + a m a- m ) , 



(28) 



the Hamiltonian can be rewritten 

H r =H [P°]+ H m [a m ,a; 



>a*l (29) 



meN, m#0 



Notice that the Hamiltonian vector field X^ on T a is just 

on To and on 
In other words, for any state in T a , 



the sum of the Hamiltonian vector fields X^ on io^u.i u 



r, with H := £ 



meN-(0| 



H„ 



the time evolution can be deduced by composing the evolution 
of the projections in To and in f , where the dynamics are dic- 
tated, respectively, by Ho and H. In particular, we see that the 
unitary implementability of the dynamics in T a at the quan- 
tum level depends only on whether the finite transformations 
generated by H on f can be unitarily implemented. 

The covariant phase space V, on the other hand, can be de- 
scribed using as coordinates the canonical pair (qo,po) and 
the set of (complex conjugate) annihilation and creation-like 
variables {(A„, A*)}„ e z-{0)- Let us denote the covariant phase 
space, in such a coordinate system, by Va- We can now 
separate the zero modes exactly as before, namely, Va can 
be viewed as the direct sum of the subspace V for which 
qo = Po = and the subspace Vo whose vectors have qo and 
po as the only nonzero components. 

In the following, we will respectively denote states in To 
and Vo by uq :- (0o>^) an d v 'o := (qo,Po)- Similarly, states 
in f and V will be denoted by \y m } := {(a m , a* m , a_,„, a*_ m )} and 
m m ] := {{A m ,A* m ,A_ m ,A*_ m % with meN- {0}. 

Since <p and are the configuration and momenta at S, of 
tp, we can express the coefficients <p n and P^" in terms of A„ 
and A*_ n for all n e Z - {0}, and in terms of qo and po for the 



zero mode, getting in this way a map M from Va to a . In 
addition, Eq. Ml\ defines a map M from F^pa to T a . Then 
the composition MM provides us with a map M : Va — > r„. 
A straightforward calculation shows that for the zero modes 



1 lnf 
uq - | q j | Vo , 



(30) 



where uq and Vo are treated as column vectors. For the rest of 
modes (meN- {0}), one obtains y m = M m {t)J[ m with 



MJf) = 



C m (t) d m (t)^ 

<(f) d*Jf) 

djf) c m (t) 

di{t) o o c'jt)) 



(31) 



and 



nm 

c m(t) = J—[H (mt) 



itHi(mt)] . 



djf) = J^[Hl{mf)-itHl{mt)\, (32) 

Here, i/i is the first-order Hankel function of the second kind 
l2lll . Note that the map M is such that M(Vo) = F and 
M(V) - f . Besides, the determinant of the linear transfor- 
mation J30l >. as well as that of M m , is equal to the unity. It 
hence follows that the M m 's are Bogoliubov transformations. 
Thus, we get a time-dependent canonical transformation from 
V A to r a . 

A generating function for this transformation (that depends 
on some appropriately chosen complete sets of compatible 
components -under Poisson brackets- both for Va and T a ) is 



r m (t) 



i 



-(po) lnt-p (po, 

= ia*_ m [c m (t)A- m + d m (t)A* m ] 
- ia m [d* m (t)A- m + c* m (t)A* n ] . 



(33) 



for m = and meN- {0}, respectively. After a straight- 
forward calculation we find that the partial derivative of this 
generating function with respect to its explicit dependence on 
the time coordinate t has the following form when expressed 
exclusively in terms of the components of the states in T a : 



1 



0,To = f t (Pt? , d t T m = H m [y m ] . (34) 



Therefore, we get 

d,H",y] = 2 5 '^ m[M ' r] 



H r [u,y\: 



(35) 



where H, is precisely the Hamiltonian ( I29l l. 

At this point of the discussion, it is worth recalling that, 
given a canonical transformation from certain symplectic vec- 
tor space Ei := {(g,-, p,)} to another one E2 := {{Qi, Pi)} which 
is determined by a generating function F that is explicitly time 
dependent [27], and assumed that the dynamics in E t is dic- 
tated by the Hamiltonian H\[q,p\, the corresponding Hamil- 
tonian in E 2 is H 2 [Q,P] = Hi[q{Q,P),p(Q,P)} + d,F[Q,P]. 
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Taking into account Eq. d35i . we then see that the dynami- 
cal evolution in F a , generated by H n arises entirely from the 
time dependence of the canonical transformation. As we have 
pointed out before, the total Hamihonian in Va is identically 
zero and there is no time evolution for the states (vo,{M m }). 
Obviously, this vanishing of the Hamiltonian applies as well to 
the restrictions to the subspaces Vo and V. In particular, while 
the states in f evolve along the integral curves of the Hamil- 
tonian vector field X^, the states in V are "frozen". Hence, 
an initial state \y m {to)} in F will evolve to the state \y m (t)} 
determined by the transformation y m (t) = U^^(t,to)y m (to), 
where 4 U^^(t, to) := M m (?)M m (?o) _1 . In contrast, the corre- 
sponding states in V, specified by M m (to) = M m (toY l y m {to) 
and M m (t) = M,„(ty l y m (t), will be related via the identity 
map, so that they actually coincide. 

In coordinates J[ m rather than y m , the finite transformation 
U^ M is given by U ( ™> m (t, f ) = M m (t y l M m (t). As a result, 
the complex antilinear part of the finite transformation gener- 
ated by H — H r — Hq in V is given (for each m e N - {0}) by 



D m (t, t ) = — [toHl(mto)H*(mt) - tW Q (mto)H\(mt)\ , (36) 

which is not square summable in m, as has been proved in 
Ref. Jg]. Therefore, the finite transformation provided by 
M m (to)~ l M m (f) (with m running in N - {0}) cannot be uni- 
tarily implemented. Moreover, since the antilinear part of 
\ 7 ](t, to) differs from that of ^(t, to) just by a sign in 
the phase of the coefficient c m (t), as one can easily check, we 
see that the finite transformation generated by H in T is not 
unitarily implementable. It is worth emphasizing that, how- 
ever, this it is not the case for the dynamics in V; indeed, since 
such a dynamics is generated by the zero Hamiltonian, the 
evolution is described by the identity transformation, which is 
of course unitary. 

Actually, since U^'^(t, to) is just a composition of the Bo- 
goliubov transformations M m {t) and M m (to), the lack of a uni- 
tary implementation of the dynamics in f follows from the fact 
that the antilinear part of M m (t) fails to be square summable 
for generic t > 0. Hence, whether or not the dynamics can 
be unitarily implemented depends entirely on the behavior of 
d,„(f) for large integers m. This depends in turn on the Hamil- 
tonian via Eq. d34i . which relates the generating function of 
M m (t) with the generator of the dynamics after the canonical 
transformation has been performed. Our analysis about the 
lack of unitarity will therefore focus on the identification of 
those characteristics of the Hamiltonian that are in the origin 
of the failure of square summability. In doing so, we will be 
able to establish how critical this problem is and whether it 
can be corrected with small modifications to the dynamics. 
Roughly speaking, we will be able to determine how far the 
considered evolution is from being unitarily implementable. 



With this aim, in the next subsection we introduce a correc- 
tion to the Hamiltonian H, which might be viewed as a per- 
turbation or a quantum modification, and discuss the square 
summability of the antilinear part of the transformation gen- 
erated by the new Hamiltonian. 



B. Modified Hamiltonian 

Motivated by our previous analysis, let us assume now that 
a certain linear (free) field theory can be described by either of 
the two symplectic vector spaces Ea := {^M m } or Eb := {S„,}, 
where W m := (A m ,A* m ,A- m ,A*_ m ), <B m := (B m , B* m , B. m ,B*_J, 
and m e N— {0}. Here, (A m ,A* n ) and (B m , B* m ) are annihilation 
and creation-like pairs. In addition, let us suppose that the 
canonical map M m (t) from Ea to Eg, which in general may 
depend on the time coordinate f, is a Bogoliubov transforma- 
tion of the form jilt , with |c m (f)l 2 _ \d m (t)\ 2 = 1. Besides, we 
assume that the total Hamiltonian in Ea is zero, so that the 
Hamiltonian in E B is given by H — 2 m£ N-{0| d t Tm, where T m 
is a generating function of the transformation M m (t). Further- 
more, we admit that (in coordinates J[ m ) the partial derivative 
of the generating function T m with respect to its explicit time 
dependence is 

d,T m = H m [$t\ : = 2m [pia(x)A m A-,„ + /j.*a*(x)A* m A*_ m ] 

+ 2mAb(x)[A* m A m +A*_ m A- m ], (37) 

where x := mt (strictly speaking we should write x m , however 
we drop out the subindex to simplify the notation). Besides, A 
and fi are a real and a complex constant 5 , respectively, and 

a(x) = ^{[ff (i)] 2 + [H,(i)] 2 }, 

b(x) = !4# (*)| 2 + |#i(*)| 2 }. (38) 

At this stage, it is convenient to point out the analogy with 
our symmetry-reduced model. The Bogoliubov transforma- 
tion M m (t) has the same form as that in Eq. J3 II . The sym- 
plectic vector spaces Ea and Eb play the role of V and T, re- 
spectively. In this sense, note that the total Hamiltonian in Ea 
is zero. Moreover, setting A = fi = 1 in equation fl37i one 
merely gets the H m contribution to the total Hamiltonian H 
in F, expressed in coordinates !A m 0]. Hence, we can think 
in terms of the spaces V and T, and regard the phase space 
function H as a modification of the Hamiltonian H. Defining 
p := A — 1 and e = fx — 1, one may view the case \p\ <k 1, 
|e| <K 1 as a perturbation of the Hamiltonian, arising from 
certain (classical or quantum) corrections to the dynamics. 

As we already know, the unitary implementability of the 
dynamics dictated by H in Eg depends on the square summa- 
bility of the antilinear part of M m (t). By analyzing this 



5 We might allow for an ^-dependence in A and fi, but this would unnecessar- 
ily complicate our discussion. Some comments about this generalization of 
4 Actually J/W ,(fo,i ) = land U^'\ ,(t 3 ,fi) = U™. :(t 3 ,t2)Uf\ fi 2 ,h). the analysis are presented in Subsec. III.C. 
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summability we will relate the failure of unitarity with the pre- 
cise form of H. 

We first determine the relations that Eq. d37l > imposes on the 
complex functions c m (t) and d m (t) that specify the Bogoliubov 
transformation M m {t). A generating function for this transfor- 
mation is 



T,„{t) = iB*_ m [c m (f)A- m + d m (t)A* m ] 
- iB m [d* m (t)A- m + c* m {t)A* m ] . 



(39) 



Taking the (explicit) time derivative and using then the inverse 
of M m {t) [which is easily calculated from Eq. d3 lil. one ar- 
rives at the following expression, exclusively in terms of the 
coordinates S m , 

d,T m = iB*_ m B- m (c m c* m - d m d* m ) + iB m B- m (c* m d* n - d* m c* m ) 
+ iB* m B m {d* n d m - c m c m ) + iB* m B*_ m (d m c 

m Cmdm)- (40) 

The dot denotes the (total) derivative with respect to the time 
coordinate t. From now on, we do not generally display the 
dependence of c m and d m on this coordinate in order to sim- 
plify the notation. 

By translating also into coordinates T> m the Hamiltonian 
H m , the condition fl37i that the dynamics in Eb arise entirely 
from the time derivative of our canonical transformation can 
be seen to reduce to the following system of first-order (com- 
plex) differential equations for c m and d m : 

= ^(c^c m -^t/ m ) + m^(x)(|c„,| 2 + |c/ m | 2 ) 

- m \p.a(x)c* m d m + ^ a{x)c m d* m ] , (41) 
= - (d m c m - c m d m ) + 2mAb(x)c m d m 

— m \pLa(x)d m d m + ^i*a"{x)c m c ln \ . (42) 

Let us call Y m the ratio d m /c m . Since \c m \ 2 - \d m \ 2 = 1, we 
have that \Y m \ is strictly smaller than the unity and, in terms of 
it, the (complex) norms of c m and d„, are 



1 



\d„ 



(43) 



1 - |F ffl | 2 ' ™ 1 - |TJ 2 

Realizing that (d„,c m - c m d,^)lc 2 m is just the time derivative of 
Y m and performing the change of variable t t-» x — mt (so 
that Y m = mdY m /dx), it is easy to see that Eq. ( 142 \ can be 
rewritten as 

i dY,„ 



2 dx 



-00 + 2Ab(x)Y m (x) - fm(x)Y^.(x) - (x*a*(x) = . (44) 



Remarkably, this differential equation for Y m is independent 
of the positive integer m (regarding x as the relevant variable). 
Using this universal character of the equation, valid for all 
values of m, we can drop out the subindex in the function Y m 
and consider it as a single function Y for all the modes of our 
system. With the convenient redefinitions 

z(x) := exp (-2ix) Y(x) , A(x) := exp(2/x) a(x) , (45) 

we then arrive at the following equation for z. 

— (x) = 2iz(x) [2Ab(x) - l]-2ifiA(x)z 2 (x)-2i^A*(x) . (46) 
dx 



In addition, given a function z satisfying Eq. J46> and re- 
membering that \c„,\ 2 - \d m \ 2 = 1, it is straightforward to see 
that the differential equation (I41> is equivalent to 



d lnc„ 

dx 



00 = 2i \jjA(x)z(x) - Ab(x)] 



(47) 



which again is a universal equation for all modes m e N - {0}. 
We suppress the subindex m and consider only one func- 
tion c(x), which can be obtained by direct integration of 
Eq. J47i [except for a multiplicative constant that can be 
fixed with an initial condition for c]. Finally, the function 
d(x) := exp (2ix)z(x)c(x) provides the missing coefficient of 
our Bogoliubov transformation, namely, d m (t) = d(x = mt). 
Using Eq. J43I and \Y(x)\ = \z(x)\, we conclude 



\d m (t)\ 2 = W(x = mt)\ 2 



\z(x = mt)\ 2 
1 - \z(x = mt)\ 2 



(48) 



Therefore, an important consequence of the observed univer- 
sality is that the square summability of the coefficients d m (t) 
at any fixed positive value of t, which is only sensitive to the 
behavior for large m, turns out to depend exclusively on the 
behavior of the function z(x) when x approaches infinity (be- 
cause x — mt grows linearly with m for all t > 0). Thus, 
to discuss the square summability of d,„, we only need to con- 
sider Eq. i46i and exploit our knowledge about the asymptotic 
behavior of the functions A(x) and b(x). 

From Hankel's asymptotic expansions of Hq and Hi l2lll 
one gets that, for x » 1, A(x) and b(x) are given by the asymp- 
totic series 6 



aoo = l - y ^l- 

4 /-J (2x) 2( - k+n 

1 - ( _iy- 
b(x) = - 2, 



ft,n+l/2 &+1/2.B+1/2 



Ax 2 



A^(2xf^ 



CTk,n + 



0£+l/2,n+l/2 



Ax 2 



(49) 



where 



frjn := (0,2fe)(0,2«)-(l,2fe)(l,2«), 

a-k, n := (0,2fe)(0,2«) + (l,2fe)(l,2n), (50) 

and (Ic, n) is the so-called Hankel symbol: 

Y[k + n+\) 



(k,ri) :-- 



n\T[k-n + \)' 



(51) 



The asymptotic series representation for b(x) contains only 
even powers in 1 jx, while the corresponding asymptotic series 
for A contains, in principle, both even and odd powers. With 



' In general, we will say that a function / admits an asymptotic se- 
ries at infinity if there exists a series of the form YZLnfkl^ such that 

lim,.^„ - 2f =0 A/ xk }\ = for a11 N > < see e -§- Ref - El)- 



9 



the change of variable y — l/x, these asymptotic expansions 
can thus be written in the form 

oo oo 

A(y) = Y j A k y k , b(y) = Y J b k y k , (52) 



k=0 



k=0 



where b%k+i - for all k e N. From Eqs. ( 1491 . we get in 
particular 



Ao = i— =0, A 1 = — = --, 



^ = 



4 ~ 2 



(53) 



Employing the asymptotic expressions for A(x) and b(x), a 
formal asymptotic series for z can be constructed. More pre- 
cisely, introducing expansions d52l in the differential equation 
J46i (with the change y = l/x) and writing z as an asymptotic 
series in y, namely z = Yik Zky k , one obtains 

oo 

zfriAo + ZbCl " 2Ab ) + ,u*A + £ y*S t = , (54) 



where, for each k > 1, 



m=0 7=0 m=0 

• a. 1 2 

+ z* + fi A k + i — — Zk-i + pA k z . 



In particular, the term independent of y in Eq. J54l > must van- 
ish. Since bo = 1/2 and An = 0, one gets (with p — A - 1) 



ZoP = 0. 



(56) 



Hence, provided that p does not vanish, the coefficient zo must 
be zero. As a consequence, the resulting function z will tend 
to zero as x — > oo, which in turn implies that \d\ vanishes in 
the limit of large values of x. 

The rest of terms in Eq. i54\ require the vanishing of S * for 
all k > I. Substituting zq =0 and the values of b(, and Ao, it is 
straightforward to derive the following recurrence relation for 
the complex coefficients of the asymptotic series of z: 



1 

Zk - - 
P 



P 



k—i m k 

V ^ A^J^ Zm-jZj) - 2(1 +p) 2^ b m Z k -„ 

m=0 7=0 m=l 



, A . .k-l i 
V A k + 1 ~^-Zk-l 



(57) 



For the first coefficient, we get z\ = /u*A*/p = -//*/(4p). 
Thus, for sufficiently large values of x, the differential equa- 
tion ( 1461 should admit a solution z such that 7 



(58) 



7 A function f(x) is o(l /je*) at infinity if lim A ^oo = 0. 



As we will see, this result suffices to prove the square summa- 
bility of the coefficients d m . Of course, the above behavior 
is not allowed for z when Ib^A - 1 = p vanishes, as it is the 
case for the Hamiltonian H. This explains the break down of 
unitarity for that specific case. 



C. Unitarity of the Modified Dynamics 

Let us finally show that the deduced asymptotic behav- 
ior for z(x), together with relation ( I48l >. guarantee the square 
summability of the sequence \d m (t)} for all fixed, strictly pos- 
itive values of the time coordinate t. From Eq. ( I58> , we see 
that, at infinity, lim MCO \xz + fi*/(4p)\ = 0. So, given any 
constant e > there exists a positive number xq(s) such that 
\xz + p*/{4p)\ < s for all x > xq(s). One can see that this 
inequality implies that, for all x > xo(e), 



2 _ R{e,p,n) 



s\p\ , \p\ 2 



IzT < ' V' ' ■ R(^,p,fj.) := 2s 2 + ^ + ^- . (59) 
x l \p\ 16p z 

Let us now choose a number xq in the interval (xo(e), oo) 
such that 770 := R(s,p,p)/x^ < 1, which is clearly always 
possible. Then, for all x > in, 



2 R(e,p,n) 



z < 



< 770 < 1 



(55) Therefore, for all x > in we have 



1 



1 - ?7o 



<4, «b 



R(s,p,p) 
1 -J7o 



(60) 



(61) 



Notice that Rq is a finite and strictly positive constant. Em- 
ploying Eq. J48I and this inequality, we obtain that for all 
m > Mq := int{io/f) + 1 (where int{x) is the integer part of x) 



\d(mt)\ 2 < R " 



(mt) 2 ' 

Therefore, we conclude that, for every fixed t > 0, 

Mo 



(62) 



meN, m#0 



m=Mo+l 



Z i«i 2 < ^i«i 2+ 5 z ^ 

m=l 

< J>(mr)| . 6f2 



: + ^ <co . (63 ) 



Thus, c/ m is square summable, and the dynamics generated by 
H in Eb is unitarily implementable. The proof, which makes 
use of Eq. J58I . fails when p — and, in particular, for the 
Hamiltonian H. We hence see that unitary implementability 
is extremely sensitive to the value of p. For instance, for every 
nonvanishing real p in any neighborhood of zero (in fact, for 
allp e R-{0}), H m [M\ = H m [J[]+2mpb(x)(A* m A m +A% n A- m ) 
gives rise to a unitarily implementable transformation. 

On the other hand, it is worth emphasizing that, to prove 
the square summability of d m , we have not actually employed 
the existence of an asymptotic series for a solution z of Eq. 
i46l . What we have used in fact is a weaker property, namely, 
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the existence of a solution with the behavior d58l >. A rigorous 
proof of this existence is given in the appendix. 

In our analysis, we have assumed a specific form for the 
Hamiltonian, in particular that A and p are a real and a com- 
plex constant. It is nonetheless possible to generalize our dis- 
cussion to other cases. Suppose, e.g., that the function Ab(x) is 
replaced by a new real function b(x) and that pa(x) is changed 
into a(x) := exp [-2/t/r(x)]A(x) (and similarly for its complex 
conjugate), with iff and A being a real and a complex function, 
respectively. Let us then call /3{x) := 2b(x) - di[r(x)/dx and 
assume that y6(x) and A(x) admit asymptotic series at infinity 
such that Ao = 0. Defining now z(x) = exp [-2ii[f(x)]Y(x), it 
is straightforward to see that the same line of reasoning pre- 
sented in the previous subsection leads to an equation analo- 
gous to Eq. i46i , but with the replacements of z by z, [2Ab- 1] 
by and pA by A. The counterpart of Eq. d56l > is then 
Zofio = 0, which implies that zo vanishes unless so does y8o, 
which plays now the role of p. In addition, the analog of Eq. 
(H7J is 



Zk = 



1 

To 
i 

To i 



k-l ( m \ k 

m=0 V, 7=0 / m= \ 



t, .k-l_ 

A,, + l —r—Zk-l 



(64) 



Therefore, we get z(x) = A*/(j3 x) + o(l/x). Again, this 
asymptotic behavior suffices to guarantee the square summa- 
bility of the antilinear part of the map M m (f) for all t > 
provided that /3o differs from zero. 

As a final comment, let us consider the formal quantum 
expression for H: 

: S :=: H : +2p ^ m b(mt) [A] n A m + P m A- m ] . (65) 

meN, mj=Q 

Because of the unitary implementability of M m (t), we know 
that d65l > generates the unitary evolution operator through 
which the basic operators B„ and B,', evolve (when these basic 
operators are represented as the annihilation and creation op- 
erators on the Hilbert space fig constructed from Eg and its 
corresponding complex structure Jb). On the other hand, even 
though the operator d65l > generates a map which acts unitarily 
on the Hilbert space Iri^ (constructed from F and its associated 
complex structure /f ), we know that this map does not corre- 
spond to the actual time evolution of the basic operators B„ 
and B\ (now represented as the annihilation and creation op- 
erators on Hi). The quantum generator comes from a phase 
space function which certainly can be considered as close as 
one wants to H (the generator of the dynamics in F), but does 
not coincide with it. Nevertheless, it is interesting to note that 
if p is regarded as a constant of quantum origin, e.g. by setting 

p proportional to Ti, then the classical limit of H would be just 
H. That is, in spite of the lack of a unitary implementation 
for H, if we consider that the modification of the Hamiltonian 
arises from a quantum correction, then we will get a unitary 
map whose generator, in the naive limit ft — > 0, provides the 
classical dynamics in F. 



IV. CONCLUSIONS AND FURTHER COMMENTS 

We have analyzed the impossibility of obtaining a unitary 
implementation of the dynamics in the polarized Gowdy T 3 
model with the quantization put forward in Ref. |5], a prob- 
lem that has recently been pointed out in Refs. |(J and 0. 
With this aim, we have first presented a complete derivation 
of the model starting with general relativity and introducing a 
symmetry-reduction and gauge-fixing procedure. Employing 
then a time-dependent map from the covariant phase space to 
the canonical phase space of the system, we have been able to 
reformulate the issue of unitary implementability of the evo- 
lution as a question about the square summability of the an- 
tilinear part of such a map. In this process, it is important to 
realize that the total Hamiltonian in the covariant phase space 
vanishes, whereas the considered map includes in an explicit 
manner all of the time variation of the system. Exploiting this 
reformulation of the unitarity problem, we have considered 
(certain types of) modifications to the dynamics and analyzed 
whether the symplectic maps associated with them are unitar- 
ily implementable. In this way, we have traced back the failure 
of unitarity to the presence of some specific contributions in 
the Hamiltonian that generates the dynamics. In addition, we 
have seen that negligibly small modifications of these contri- 
butions suffice to restore unitarity. In the rest of the section, 
we present some comments about the main results of the work. 

In our analysis, two facts have played a particularly rel- 
evant role. Firstly, as we have noticed, there is some kind 
of universality in the behavior of the Bogoliubov coefficients. 
This has allowed us to consider just one equation [namely Eq. 
J46H in order to examine the square summability of these co- 
efficients, rather than investigating an infinite number of dif- 
ferential equations, one for each mode. Secondly, to know 
whether the modified Hamiltonian is unitarily implementable, 
instead of solving the universal equation d46l >. it actually suf- 
fices to study the leading term of the function z in the asymp- 
totic limit of large values of its argument. In this sense, one 
does not need to explicitly integrate the dynamical equations. 

On the other hand, we note that solving Eq. i42l amounts 
to "diagonalizing" the total Hamiltonian by means of a time- 
dependent canonical transformation, namely, to requiring that 
the terms proportional to B m B-,„ and B* lt B*_ m vanish in the 
phase space function H 2 [S] := H[S] - d,T[S]. Had we ig- 
nored the term containing the time derivative of the Bogoli- 
ubov coefficients, we would have obtained from Eq. ( 1421 an 
algebraic quadratic equation for the ratio Y(x) of the coeffi- 
cients that leads to an instantaneous diagonalization (i.e., at 
a fixed instant of time) of the Hamiltonian H. In fact, for 
A — p = 1 so that H reduces to H, one can see that using 
this algebraic equation and the relation |c| 2 = 1 + |ii| 2 , it is 
possible to recover the instantaneous diagonalization given 
in Ref. @]. In addition, we emphasize that diagonalizing 
the Hamiltonian is equivalent to the resolution of the dynam- 
ics. Indeed, as we have seen, if we solve the universal equa- 
tion d46l >. which is equivalent to the diagonalization condition 
d42l >. then c(x) can be found by simple integration of the first- 
order differential equation ( I47L whereas d(x) is determined as 
d(x) = exp (2ix) z(x)c(x). 
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Our discussion can be extended to Hamiltonians for which 
the functions p.a(x) and Ab(x) in Eq. fl37i are replaced by more 
general functions d(x) and h(x). We have seen that this is the 
case at least if a(x) is of the form exp [-2/i/'(x)] A(x) with i[r(x) 
real, and A(x) and (the real function) ]3{x) : = 2b(x) - dift(x) / dx 
admit asymptotic series with a vanishing coefficient Ao. More 
precisely, we have proved that the dynamics generated by 
those Hamiltonians can be implemented as a unitary trans- 
formation as far as the coefficient p\> differs from zero. In this 
sense, our study provides a general treatment for Hamiltoni- 
ans of the form (I37L quadratic in the coordinates J[ m . 

We have proved that the obstructions to unitarity are ex- 
tremely sensitive to perturbations of the Hamiltonian. In the 
quantization performed in Ref. @|, one can in fact approx- 
imate the dynamics as much as desired by means of unitary 
transformations. Thus, although a thorough study of the influ- 
ence of the choice of deparametrization and complex structure 
must be performed in order to elucidate if the failure of uni- 
tarity is a phenomenon inherent to the polarized Gowdy T 3 
model, our result suggests that it may be actually possible to 
restore unitarity by considering a different choice. 

Finally, employing Eq. (1381 and the asymptotic expansions 
of the Hankel functions, it is not difficult to check that, for 
large values of f, the m-th contribution to the total Hamilto- 
nian of the polarized Gowdy T 3 model [see Eq. d37t with 
A = 1] becomes H m [Jl] m m{A* m A m + A*_ m A- m ). So, the 
Hamiltonian for asymptotically large values of t is given by 
H[M] ~ ZneZ-{0} M ^nA„, which is a combination of har- 
monic oscillators. Therefore, we could approximate the dy- 
namics in an unitary way by means of some modified Hamil- 
tonian H and, for large values of t, by a sum of harmonic 
oscillators. Actually, corrections to the dynamics are not un- 
expected, e.g., from a quantization in the framework of loop 
quantum gravity. It might happen that some kind of (fictitious) 
effective dynamics could take place as a result of the smearing 
of the initial singularity by quantum effects (like it occurs, in 
fact, for a flat isotropic universe, where the cut-off for curva- 
tures provides us with an effective Friedmann equation ll29lo . 
This effective dynamics could correspond to a perturbation of 
the evolution generated by H which might be unitarily imple- 
mentable. 
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APPENDIX A: PROOF OF THE ASYMPTOTIC BEHAVIOR 

In this appendix, we want to prove that Eq. ( 146 \ admits one 
solution which, at infinity, has the asymptotic behavior (j58jl, 
provided that p + 0. Let us start by defining a new function 



w(x) by means of the relation 



z(x) 



w(x) 



4px 



(Al) 



Substituting this expression in Eq. ( I46l l we obtain an equiva- 
lent nonlinear differential equation of the Riccati type 



dw A(x) 2 

— (x) = w(x)/3(x) - 2ifx w (x) + a(x) , 

dx x 



(A2) 



where 



fj(x) :- 2i [2Ab{x) - 1] H v i , 

x p x 

a(x) :- — 2ifi*x[A*(x) — A\] — i- — [b(x) - bo] 



It . At*M 2 A(jc) 
Apx 8p 2 x 



(A3) 



The constants A] and bo are given in Eq. d53l > and we have 
used A = 1 + p. 

In order to arrive at the desired result about the asymptotic 
behavior of z(x), we only have to demonstrate that Eq. iA2l 
admits a solution that tends to zero at infinity. 

Employing the asymptotic expansions d52t of the functions 
A(x) and b(x), recalling that Ao = b\ — 0, and making use of 
Eq. fl49i to compute the coefficient A2 = i/16, one can rewrite 



1 - 

P(x) = 2ip+-+ fi(x) , 
x 



a(x) 



1 l\p* 



2 pj Ax 



+ a(x) , 



(A4) 



with a(x) and y6(x) being (9(1 /x 2 ) at infinity [we say that a 
function fix) is (9(1 /x") at infinity if \x"f(x)\ admits a finite 
limit when x — > 00]. Explicitly, these functions are 



fj(x) := AiA [b(x) - b ] + i 
a(x) : 



M 2 A(x) 



:= —2ip*x 



p x 

K 



A*(x) - A* - 



- i — [b(x) - bo] 
P 



^*M 2 A(j) 
8p 2 x 



(A5) 



On the other hand, the function A(x)/x that multiplies w 2 (x) 
in Eq. flA2b is also (9(1 /x 2 ) asymptotically. For solutions w(x) 
that are small at infinity, we then expect the quadratic term 
in our Riccati equation to be negligible. We will hence ap- 
proximate our equation by a linear one that can be explicitly 
solved, find for it a solution that tends to zero at infinity, and 
prove that, for that solution, the removed quadratic term can 
in fact be neglected in the original differential equation. 
For the linear differential equation 



dwi 

—(x) = Wi(x)J3(x) + a(x) 
dx 



(A6) 
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all solutions can be constructed starting with those of the as- 
sociated homogeneous equation. Using Eq. iA4l . the homo- 
geneous solutions can be found to be proportional to 



wy(x) = xexp 



•J \X 



2ipx + | dx /3(x) 



(A7) 



Note that the asymptotic behavior of /3 guarantees that the in- 
tegral that appears in this expression is well-defined. Solu- 
tions to Eq. ( IA6> [modulo the possible addition of w'j(x) times 
a complex constant] are then of the form 



w/(x|x ) = w, (x) 



f x ,_a{x) 
I ax 

Jxc vvf(x) 



(A8) 



where x is a constant. A convenient integration by parts leads 
then to 



wi(x\x ) 



J. l 



+ <(x) 



dx 



vt^(x) dx 
a(x) 



a(x) 



{2ip+0(x)} 



rf{x){2ip+P(xj\ 



(A9) 



Remembering that a(x), /3(x), and wj(x)/x are respectively 
0(1 /x), 0(1 /x 2 ) and 0(1) at infinity, it is possible to see that 



the integrand in the above expression is 0(1 /x 5 ). Therefore, 
the integral converges when xq tends to infinity. 
In that limit, one gets the particular solution 



w l (x) - w\ (x) I dx 



f 



w'j(x) dx 



a(x) 



[lip + fax)} 



a(x) 



lip + P(x) 



(A10) 



One can see (e.g. using L'Hopital's rule for the term contain- 
ing the integral) that this solution tends to zero when x — > oo. 
Furthermore, repeating the explained procedure of integration 
by parts, one can show that the total contribution to wf (x) 
coming from the factor that includes the integral is 0(1 /x 2 ). 
Using Eq. JA4> . one then concludes that the asymptotic be- 
havior of wf (x) is 



wf(x) 



// 11 

ho - 

8/px \x 



(All) 



It is now a simple exercise to check that, in the Riccati equa- 
tion JA2t . the quadratic term is o(l/x 3 ) at infinity for the so- 
lution wf (x), which is in fact negligible when compared with 
the rest of terms in the equation (in particular with dw/dx), 
which are at least of order 1/x 2 . This concludes our proof. 
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